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We study the effect of modifications to General Relativity on large scale weak lensing observables. 
In particular, we consider three modified gravity scenarios: f{R) gravity, the DGP model, and TeVeS 
theory. Weak lensing is sensitive to the growth of structure and the relation between matter and 
gravitational potentials, both of which will in general be affected by modified gravity. Restricting 
ourselves to linear scales, we compare the predictions for galaxy-shear and shear-shear correlations 
of each modified gravity cosmology to those of an effective Dark Energy cosmology with the same 
expansion history. In this way, the effects of modified gravity on the growth of perturbations are 

00 . separated from the expansion history. We also propose a test which isolates the matter-potential 
' relation from the growth factor and matter power spectrum. For all three modified gravity models, 
I the predictions for galaxy and shear correlations will be discernible from those of Dark Energy 

OA . with very high significance in future weak lensing surveys. Furthermore, each model predicts a 

measurably distinct scale dependence and redshift evolution of galaxy and shear correlations, which 
can be traced back to the physical foundations of each model. We show that the signal-to-noise for 
detecting signatures of modified gravity is much higher for weak lensing observables as compared to 
the ISW effect, measured via the galaxy-CMB cross-correlation. 
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A number of independent observations, ranging from Supernovae to the cosmic microwave background (CMB) and 
, its cross-correlation with foreground galaxies [1-5], have now firmly established that the expansion of the Universe 
J> ■ is accelerating, and that a purely matter-dominated universe is not consistent with measurements. Commonly, these 
04 ' observations are ascribed to an additional smooth stress-energy component, "Dark Energy" (DE), pervading the 
Universe [6]. However, instead of attributing the accelerated expansion to our lack of understanding of the constituents 
of the Universe, one can take the alternative approach and attribute it to our lack of understanding of gravity on 
cosmological scales. For example, an effective weakening of gravity on the largest scales could explain the accelerated 
' expansion and other cosmological observations without invoking an additional smooth "dark" component. Several 
such modified gravity scenarios have been proposed. The difficulty in constructing a consistent theory of gravity which 
preserves the success of General Relativity (GR) in the Solar System and the early Universe has however limited the 
number of proposed theories to only a handful. 

Seen from a fundamental physics point of view, the two approaches of General Relativity coupled with a smooth 
Dark Energy (GR-I-DE) and modified gravity are apparently quite distinct. However, there is enough freedom in both 
scenarios to match any given expansion history of the Universe. Thus, the expansion history, measured primarily with 
Supernovae and the CMB, is not sufficient to distinguish between these two "fundamentally different" scenarios. 

Going beyond the smooth background Universe offers a multitude of additional observables probing the evolution of 
structure formation in the Universe. For the purpose of constraining modified gravity, the evolution of the cosmological 
gravitational potentials, and their relation to the matter overdensities and velocities are the most sensitive probes [7- 
10]. Several such observational tests of gravity have been proposed: the Integrated Sachs- Wolfe (ISW) effect [11-15] 
directly probes the evolution of the potentials on horizon-size scales. While quite sensitive as a probe of gravity, it is 
restricted to large scales, and offers only a limited amount of signal-to-noise. The matter power spectrum probes the 
growth of matter perturbations on a range of scales. While it can be measured to high precision and to very small 
scales, and hence offers a large amount of information, it is affected by the non-linearities of structure formation, and 
galaxy bias, which are not well understood for modified gravity models. A cleaner method of probing gravity with 
galaxies is using velocity correlations [7, 16, 17]. In principle, these tests can isolate signatures of modified gravity. 
However, they will be challenging to perform observationally. 
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Several deep, wide-field galaxy surveys are in the planning stage which are designed to measure galaxy counts 
and weak lensing shear with unprecedented accuracy. Galaxy and shear correlations can be used to measure the 
expansion history and growth of structure in the Universe to high precision. Additionally, the potential of weak 
lensing observables to probe gravity has been shown in [7, 18-21]. In this paper, we investigate the constraints on 
modified gravity theories to be expected from galaxy and shear correlations in future surveys. Since weak lensing is 
caused by the potential wells along the line of sight, it can be used to measure the scale dependence and cosmological 
evolution of these potentials. In addition, correlating foreground galaxies with the shear of background galaxies offers 
a test of the matter-potential relation. We will show that future weak lensing surveys will place very tight constraints 
on gravity. In other words, they will perform stringent tests of the smooth Dark Energy scenario, which will have to 
match all of the weak lensing observables in order to remain a viable model. 

Regarding the potential of these precision tests, it is important to keep several caveats in mind: first, the non-linear 
evolution in modified gravity models is not yet understood in a realistic cosmological context. Hence, in this paper we 
will limit ourselves to scales where linear theory is appropriate. Second, correlations involving galaxy number counts 
are affected by the galaxy bias, which is a priori unknown. We will point out how the effects of bias can hopefully 
be disentangled from modified gravity effects. Third, it is always possible to find a general, non-smooth Dark Energy 
model which mimics the predictions of any given modified gravity model [22, 23]; note that such a DE model will 
have significant density and anisotropic stress perturbations on sub-horizon scales. 

In this paper, we will study three different modified gravity scenarios, focusing on their predictions for lensing 
observables: (i) /(i?) gravity [24-28]; (ii) the Dvali-Gabadadze-Porrati (DGP) model [29-31]; and (iii) Tensor- Vector- 
Scalar theory (TcVcS) [32]. While the first two models arc able to achieve accelerated expansion without dark energy, 
TeVeS is designed to explain observations without dark matter by showing MOND-behavior in certain regimes. These 
three theories encompass very different approaches to the problem of constructing a consistent theory of gravity. 
Hence, it is interesting to determine not only whether observations can detect departures from GR, but also whether 
they would be able to discriminate among different modified gravity scenarios. 

In addition to predicting a late-time accelerated expansion of the Universe, all these models have been designed to 
approach GR in the high curvature regime, which applies to the Solar System as well as the early Universe. In this 
way, the models pass Solar System tests while making predictions close to the standard cosmological model for the 
CMB and Big Bang nucleosynthesis. 

Weak lensing probes the growth of large-scale structure which is sensitive to the background expansion history, as 
the expansion rate enters the evolution equations of perturbations. Since two of the models discussed here (DGP 
and TcVcS) predict expansion histories somewhat different from ACDM, wc compare them with smooth Dark Energy 
models designed to mimic the expansion history of these models by employing a suitable equation of state w{a). 
Deviations in the weak lensing predictions of DGP and TeVeS from the corresponding dark energy model are then 
solely due to modified gravity effects on the growth of structure and the matter-potential relation, and independent 
of the exact expansion history. 

The structure of the paper is as follows. In Section II, general expressions for galaxy and shear correlations applicable 
to modified gravity are given, and we point out out the different channels through which modified gravity can affect the 
observables. We also present a cosmological test of the Poisson equation using weak lensing and galaxy correlations. 
In section III, we briefly describe the modified gravity models studied here, and outline the characteristics which 
determine their lensing predictions. We then present forecasts for galaxy-shear and shear-shear correlations in future 
surveys in section IV, showing the deviations of the modified gravity predictions from those of GR-I-DE, as well as 
the expected signal-to-noise in actual surveys. We conclude in section V. 

II. WEAK LENSING CORRELATIONS 

A. Weak lensing in modified gravity 

In this section, we summarize the lensing observables considered in this paper, generalizing the standard expressions 
to the case of modifled gravity. We assume a flat Universe throughout. Our metric convention is [33]: 

ds'^ = a'^irf) [-(1 + 2'^)dTf + (1 + 2$)dx2] , (1) 

where 0(77) is the scale factor, -q denotes conformal time, and \E' are the cosmological potentials. In General 
Relativity (GR), $ = in the absence of anisotropic stress, which is the case at late times in ACDM and most 
dark energy models. In addition, the potentials arc conventionally re-cxprcssed in terms of the matter ovcrdcnsity 
using the Poisson equation. In case of modified gravity, both of these assumptions do not necessarily hold. In the 
following we give generalized expressions for the different observables which apply to the case of modifled gravity (see 
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also [7, 18]). We restrict ourselves to linear perturbations valid on large scales (k < 0.1 h/Mpc). Since the f{R) 
and DGP models recover the GR+DE limit at early times in the matter-dominated epochs, it is convenient to refer 
observables to the matter overdensity at high redshift Zm- For definiteness, we take Zm = 50. Then, for any viable 
f{R) and DGP model, P{k, Zm) is identical to that expected in GR. 

Weak lensing surveys use the observed ellipticities of galaxies to reconstruct a map of the cosmic shear, which 
can then be used to infer the convergence k [34, 35]. In the following, we will use the terms convergence and 
shear interchangeably, understanding that they are reconstructed from observed galaxy ellipticities. Using standard 
approximations, the convergence in the direction n, given by the Laplacian of the projected lensing potential, is 
expressed as a line of sight integral to the source at redshift Zg'. 

where x(z) is the comoving distance out to redshift z, H(z) = ^/a is the expansion rate, = l/2($ — \E'), 
denotes the Laplacian in terms of comoving coordinates, and the lensing weight function is given by: 

WLiXs,x) = —iXs-x) (3) 

Xs 

Hence, the convergence is sensitive to the potential in a range of distances centered around x ~ Xs/2- Note that 
in the case of GR and smooth Dark Energy (referred to as "GR+DE" in the following), = $ = — ^f, and we can 
apply the Poisson equation: 

47rGaV„ 6 = ^QmH^a-'d (4) 

In modified gravity, this equation does not hold anymore, and we define a generalized "Poisson factor" relating 
the potential to the matter overdensity S in Fourier space: 

Apparently, a departure of the Poisson factor D^_ from its expected value in General Relativity is a signature of 
the modification of the Poisson equation in alternative gravity theories (for a discussion of general, non-smooth Dark 
Energy models, see section III D). In addition, the growth of potential and matter perturbations is affected by modified 

gravity. The growth factor of matter, which we normalize at z = z„i = 50 so that Dm{k, z) = S{k, z)/5{k, z = z„i), 
becomes scale-dependent in many modified gravity models. We can combine the growth of structure and the Poisson 
equation to obtain: 

^-^'''^^ -D^_(k,z)Drr,{k,z). (6) 



5(fc,^; = Zm) 



This allows us to express the galaxy-galaxy, galaxy-shear and shear-shear correlation coefficients, C^^{t), C^'^{£), 
C""*(£) as follows: 

C«*(f) = dkk^P{k,Zm) j dzbiWg,{z)Dm{k,z)Mkx{z)) j dz'bjWg^{z')Dm{k,z')Mkx{z')) (7) 

= I Jdkk^P{k,Zm) J dzbiWg^{z)Dm{k,z)je{kx{z)) J dz'W,^{z')k^D^_{k,z')Dm{k,z')u{kx{z')) (8) 

C-i-i^i) = dkk^P{k,Zm) J dzW^,{z)k''D^_{k,z)Dm{k,z)je{kx{z)) 

X J dz' W^^ {z')k'D^_ (k, z')Dm{k, z')jeikxiz')) (9) 

Here, P{k, z„i) is the matter power spectrum at early times z = z„i, which in the case of the /(i?) and DGP models is 
the same as for a GR + smooth DE model with the same expansion history. The indices i, j denote different redshift 
bins, i.e., Wg-{z) is the galaxy redshift distribution for bin i, normalized to 1, hi is the galaxy bias for the same bin, 
and the shear weighting function of redshift bin j, W^j, is given by: 

W., {z) =j^J^^ dzMxizs), x{z))Wg^ {zs) (10) 
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For I > 10, we can apply the Limber approximation [36], which simpUfies the expressions considerably: 

C^^^^ii) = ldzW.hW„{z)b,W,,{z)[Dl{k,z)P{k,Zm)]^^i^ (11) 

Co^-^ii) = j dz^hWg,{z)W,,{z)[D^_{k,z)Dl{k,z)eP{k,Zrn)]k^i+m (12) 

(£) = / dz^^W,^ {Z)W,^ (Z) \dI_ (fc, Z) Dl^ik, Z) k'P{k, Zra)] ,,,,, (13) 



. 1 + 1/2 



Weak lensing correlations thus depend on cosmology through three distinct channels: (i) the background expansion 
history, via H{z) and Wl{Xs,x{z))^ (ii) the growth of perturbations and the Poisson equation, through Dm{k,z) and 
D<i,_{k,z); and (iii) the matter power spectrum at early times, P{k,Zm)- In case General Relativity holds and the 
accelerated expansion is due a smooth Dark Energy, only the expansion history (i) is affected by the Dark Energy, i.e. 
modified from ACDM. In case of the late-time acceleration modified gravity models f{R) and DGP, both channels 
(i) and (ii) are affected. And the TeVeS model without dark matter modifies (i), (ii), as well as (iii) from the ACDM 
case. 

In this paper, we want to illuminate how weak lensing correlations arc sensitive to these three channels, both in 
redshift evolution and in scale dependence. Clearly, it is crucial to separate these effects in order to distinguish 
between GR+DE and modified gravity; in addition, testing Dm{k,z), D<^_{k,z), and P{k,Zm) separately will allow 
to distinguish between different modified gravity scenarios. 

The galaxy-shear correlation is proportional to the a priori unknown bias of the foreground galaxies. Assuming 
that the bias is scale-independent on linear scales, one can marginalize over this parameter for a given galaxy sample. 
An alternative is to consider the reduced galaxy-shear correlation W^{t), defined by: 

TiSi'^i m = ± (14) 

For redshift bins that are not too wide, the linear galaxy bias drops out of this expression, so that i?^" is independent 
of hi. Since it is not straightforward to compute the errors on this quantity, we will for simplicity only consider the 
galaxy-shear correlation divided by the bias, Ca'^{t)/h, in this paper. In any case, we do not expect the constraints 
to be degraded very significantly when considering instead of C^*^. 

On small scales fc ~ 0.1 /i/Mpc and above, the density contrast grows to order unity, and non-linearities in the 
formation of structure become important. Furthermore, all of the modified gravity theories studied here employ non- 
linear equations of motion which are expected to restore gravity to General Relativity in high-density environments 
[32, 37-39]. These mechanisms have been studied in special cases, e.g. spherically symmetric and/or static cases 
[32, 37]; however, the details of cosmological structure formation in the non- linear regime have not been studied in 
the case of modified gravity. Hence, we will restrict our discussion to linear scales in this paper. 



B. A cosmological probe of the Poisson equation 

Recently, techniques have been proposed to isolate the geometrical factors in weak lensing correlations from the 
growth of structure, bias, and matter power spectrum (shear ratios, or "cosmography" [40, 41]). Because of the 
separation of geometry from growth effects, cosmographic techniques apply equally well in modified gravity theories, 
i.e. when using the generahzed expressions (11)-(13). 

A similar approach as for cosmography can be used to separate out the effect of a modified Poisson equation from 
the growth factor and matter power spectrum. Observing a modification of the Poisson equation on cosmological 
scales would constitute a "smoking gun" for deviations from the GR -|- smooth Dark Energy scenario. Consider a 
very narrow galaxy redshift distribution centered around Zf. Then, the galaxy-galaxy and galaxy-shear correlation 
coefficients become: 

C«s(^) ^ E^b'^[Dl{k,z)P{k)\^^^^, Xf^xizf) (15) 
Xf 

Ca-{i) B^bW,^,^) [D^_ {k, z) D^k, z) k^P{k)\ (16) 

Xf X, 

With this simplifying assumption, ratios of weak lensing correlations can be used to isolate the effect of a modified 
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Poisson equation: 

These ratios divide out the dependency on the matter power spectrum and growth factor; they still depend on the 

expansion history, although not very strongly, via xj the Icnsing weight function W,i{zj). On linear scales, T'{t) 
is independent of scale in GR+DE, due to the scale-independence of both k^D^_ and bias b. In contrast, a scale- 
dependence is expected for gravity theories which modify the Poisson equation. Hence, measuring V{i) as a function 
of scale £ would constitute a robust test of the Poisson equation on large scales. See section III D for a discussion of 
the effect of non-smooth Dark Energy. 



III. MODIFIED GRAVITY MODELS 



In this section, we briefly present the modified gravity models and parameters adopted in this paper, and discuss 
the qualitative expectations for their lensing predictions. First, let us discuss general constraints applicable to any 
theory of gravity. In any metric theory of gravity that conserves energy-momentum, the evolution of super-horizon 
perturbations is essentially defined by the expansion history alone [10]. This can be understood in the "separate 
universes" picture: a super-horizon curvature perturbation behaves essentially as a separate Priedmann-Robertson- 
Walker Universe with a constant small curvature. Hence, the evohition of super-horizon curvature perturbations in 
modified gravity models is the same as in GR (given the same expansion history) [22, 23]. Once a scale-independent 
super-horizon relation between the metric potentials $ and ^' has been supplied by a theory, the evolution of \E' 
is fixed on super-horizon scales. In the opposite limit, when a mode is well within the horizon, one can apply the 
quasi-static approximation, neglecting time derivatives with respect to spatial gradients. In this regime, the potentials 
$, \l/ are given by a modified Poisson equation which again is determined by the theory. 

Recently, Hu and Sawicki [22] presented a parametrization of modified gravity theories which is based on the super- 
horizon and quasi-static regimes, with a suitable interpolation between both. They showed that this "Parametrized 
Post-Friedmann" (PPF) approach reproduces the predictions of the late-time acceleration f{R) and DGP models 
well, with only few model-dependent functions. Several alternative parametrizations of modified gravity have been 
proposed [7-9, 17, 23]. The PPF parametrization is most useful for us in that it describes the f{R) and DGP models 
with suitable parameters. 

The two main model-dependent quantities with impact on cosmological observables are (in the notation of [22]), 
(i) the metric ratio g{k,z), given by: 



$ -)- \|/ 

g{k,z) 



(18) 

k.z 



and (ii) the rescaling fci^) of the Newton constant in the Poisson equation [equation (4)] on sub-horizon scales: 

1 + jg{z) 

Note that for GR coupled with a smooth Dark Energy, g ~ ^ at late times. Thus, fc influences directly, 
by rescaling the GR value by 1/(1 + fa). The metric ratio g influences the evolution of cosmological perturbations 
[22, 23]: if 3 > 0, i.e. $ -|- \E' > 0, the potentials decay faster during the epoch of onsetting acceleration than in the 
GR limit. Conversely, for gr < the potential decay is slowed down and, for sufficiently negative g, inverted into a 
potential growth. From equations (11) (13), we expect this suppressed or enhanced growth of perturbations to lead 
to potentially observable effects in the weak lensing correlations. 



A. f{R) model 

In f{R) gravity, the Einstein-Hilbert action is modified by adding a certain function of the Ricci scalar, f{R), to 
the gravitational part of the Lagrangian [24-28, 42 45] : 



L_Jd''x^[R + m] (20) 
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FIG. 1: ie/i panel: Deviation of the growth factor Dm{k, z) from a GR+DE model with the same expansion history, for three 
different modified gravity models: f{R) (red/dotted; using a value of Bo = 0.4), DGP (blue/short-dashed), TeVeS (green/long- 
dashed; normalized to 1 at 2 = 5 for clarity). The three lines show different scales: k = 10~^ /i/Mpc (thick), k — 0.01 ft/Mpc 
(medium), and k = 0.1 /i/Mpc (thin). Right panel: The Poisson factor D$_ [equation (5)] scaled by k^ /Hq/^I + z), which 
reduces to 3Slm/2 in the case of an unmodified Poisson equation, as a function of redshift for ACDM and modified gravity 
models. The lines correspond to the same scales as in the left panel. 



In the case of f{R) = C = const., we recover ACDM, i.e. GR vifith a cosmological constant A = — C/2. A linear 
function ,f{R), so that //? = df /dR is a constant, simply amounts to a rescaling of Newton's constant, G ^ Gj + 
Thus, a nontrivial modification of gravity is linked to a nonzero second derivative of /. The dimensionless field 
//{ = df /dR then appears as an additional scalar degree of freedom in the modified Einstein equations. 

Different functional forms of f{R) have been proposed in the literature. They generally satisfy the asymptotic 
relation f{R) for R — > oo, so that GR is in principle restored in the high-curvature limit, which applies to the 
Solar System as well as the early Universe. Note that Solar System constraints strongly limit the possible choices of 
/ [38, 42, 45]. In addition, d^f/dR^ should be positive in order to achieve a stable high-curvature limit [26, 43]. For 
this paper, we adopt the implicit definition of f{R) proposed in [26]. Since the background curvature R is fixed as a 
function of a for a given expansion history a{t), the function /(i?(a)) can be determined from the modified Friedmann 
equations, which are obtained by applying the modified Einstein equations to the FRW metric. We assume the 
expansion history of our fiducial AGDM model, defined by fl„i = 0.27, Ha = 0.73, fib = 0.046, h = 0.7, Ug — 0.95. For 
the power spectrum normalization, we use the CMB normalization of the primordial curvature perturbations given 
by i5<^ = 4.58 • 10~^ at fc = 0.05 Mpc~^, as in the case of the ACDM model. This is motivated by the fact that this 
f{R) model makes definite predictions for the CMB (i.e., identical to ACDM), whereas measurements of the power 
spectrum amplitude today, e.g. using the power spectrum of galaxies, are affected by gravitational non-linearities, a 
regime where definite predictions for modified gravity are still outstanding. 

The amplitude of the field fpi, and correspondingly of the modification to gravity, is most conveniently parametrized 
in terms of the parameter Bq — B{a — I), where B{a) is a fimction proportional to the second derivative of /: 



B(a) 



dy/dR^ H 

1 + fR H' 



(21) 



where R{a) is the Ricci scalar of the FRW background, which in turn can be expressed in terms of the Hubble rate 
H{a) and its derivatives, and primes denote derivatives with respect to In a. Thus, the evolution of f{R) in the 
background is fixed by the expansion history together with Bq, and we then proceed to calculate the evolution of 
matter perturbations and potentials using the PPF parametrization presented in [22]. 

For a positive value of B, necessary for stability in the high-curvature limit, f{R) models generically predict a 
negative $ -|- 4", i.e., a metric ratio g < 0. On super-horizon scales, g is of order —B, and grows to a value of —1/3 on 
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small scales. The negative potential ratio leads to a slower decay of the potentials during the acceleration epoch, and 
hence a smaller growth suppression of matter perturbations (Fig. 1, left panel) with respect to ACDM. Having fixed 
the fluctuation amplitude at early times through the CMB, we thus expect an increased weak lensing signal from the 
f{R) model. 

As pointed out above, the gravitational constant is rescaled by 1/(1 + fn) in f{R) models. Since f{R) < to 
achieve acceleration, Gcff > G in f{R) gravity, and hence is larger than its ACDM value (right panel of Fig. 1). 
However, this effect is quite small, at the percent level, for observationally allowed models. For the weak lensing 
forecasts below, we adopt a value of Bq = 0.1, leading to a field amplitude today of fn{a = 1) = —0.017. For purposes 
of presentation, we use the larger value of Bq = 0.4 for the plots in Fig. 1. 

B. DGP model 

In the Dvali-Gabadadze-Porrati model [29-31], all matter and non-gravitational interactions are confined to a 

(3+l)-dimensional "brane" embedded in (4+l)-dimensional space. Only gravity is five-dimensional, and the effects 
of the large extra dimension become noticeable at the crossover scale Tc = G'^-* /2G''^\ where G^^^ is the gravitational 
constant in five dimensions, and G^^' is the usual four-dimensional gravitational constant. The crossover scale rc is 
the only free parameter of the model. One of the two possible branches of the model naturally leads to an accelerated 
expansion when the horizon scale c/H{z) reaches Vc- The expansion history in the DGP model is somewhat different 
than ACDM. It is given by [30]: 

H{z) = Ho (^V^c + V^rc + n^a-^^ , (22) 

where J7rc = l/(4?'^i?o), and, for a flat Universe, fire = (1 — f^m)^/4. Here, we adopt flm = 0.25 and h = 0.66, 
as determined from fitting Supernova distances and the CMB [14], corresponding to TcHq « 1.78. For the same 
reasons as in the case of the f{R) model, we assume a primordial power spectrum normalization of (5^ = 4.58 • 10~^ 
at fc = 0.05 Mpc~^ for the DGP model. For the evolution of perturbations in DGP, effects from the 5D nature of 
gravity, namely perturbations in the extrinsic curvature of the brane, have to be taken into account in order to close 
the perturbation equations [46, 47]. Here, we again employ the interpolation based on the PPF parametrization given 
in [22]. In order to compare the DGP predictions with those from GR-I-DE, we use a Dark Energy model with an 
equation of state Wca{a) which exactly mimics the expansion history of the DGP model. 

The DGP model has a positive metric ratio g which grows proportional to 1/{H rc) and is of order unity at z = 0; 
it is only weakly scale-dependent for sub- horizon scales. The positive g leads to an enhanced potential decay with 
respect to GR-I-DE [14], as is apparent in Fig. 1 (left panel). Correspondingly, we expect a reduced lensing signal 
for the DGP model. Note that the decay of the potentials sets in at considerably higher redshifts, z ^ 5, than in 
GR-I-DE. This is because w^b ~ —1/2 in the DGP model for z > 1 [48]. The Poisson equation is not modified in the 
DGP model. However, as the Poisson ratio 7'(£) is proportional to flm [equation (17)], the prediction for the DGP 
model shows a scale-independent departure due to the slightly smaller value of adopted for this model (right panel 
of Fig. 1). In contrast to the simple DGP model considered here, generalized braneworld-inspired modified gravity 
models can show a modified Poisson equation [49] . 

C. TeVeS 

The TeVeS model [32] is a relativistic metric theory of gravity which, as the models presented above, reduces to 
General Relativity in the high-density, high acceleration regime. However, in the high-density, weak acceleration 
regime, it behaves like MOdified Newtonian Dynamics (MOND, [50]). The MOND force law can explain galactic 
rotation curves without invoking dark matter. TeVeS generalizes this idea by adding an additional vector field (Ve) 
and scalar field (S) in such a way that the same effect of dark matter apparently seen in dynamics is also seen in 
gravitational lensing. Hence, TeVeS can be taken as an attempt to explain all cosmological observables without any 
dark matter of unknown nature. 

For this paper, we adopt the "neutrino model" from [15]. The matter content in this model is given by baryons, 
fib = 0.05, and neutrinos, ft,, = 0.17, with a cosmological constant given by flA = 0.78, in order to match the 
observed expansion history. This model has been shown to be in acceptable agrecinicint with CMB and matter power 
spectrum observations [51]. While the current TeVeS model uses a cosmological constant to achieve acceleration, 
related theories invoking vector fields have been proposed which lead to acceleration without a cosmological constant 
[52]. The evolution equations of linear cosmological perturbations in TeVeS have been derived in [51, 53-55]. In order 
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to satisfy Big Bang nucleosynthesis bounds, the TcVcS scalar field is constrained to small values throughout cosmo- 
logical history, and neither significantly contributes to the background expansion nor the growth of perturbations. In 
contrast, perturbations of the time-like vector field exhibit a growing mode which in turn boosts the growth of matter 
perturbations [54] . Hence, the perturbations of the vector field in TeVeS play the role of seeds of structure formation 
in the absence of dark matter. The vector degree of freedom is formed by the spatial components of the time-like 
vector field, which couple to the spatial part of the metric. Hence, the growing vector field only affects the potential 
which leads to a nonzero $ -|- \E' < 0, i.e. g < [54]. As in the case of the f{R) model, this produces a slower 
decay, or even growth, of the potentials at late times. However, the power spectrum at earlier times is lower than 
that of ACDM due to the considerably different history of structure formation in TeVeS. 

In addition to their effect on the potential ratio, the vector perturbations also appear on the right-hand side of the 
Poisson equation [15], and even dominate the matter perturbations on small scales. This effect can be seen in the 
"Poisson factor" (Fig. 1, right panel): on large scales (thick line), the value is significantly smaller than in ACDM 
due to the smaller Qm in TeVeS without dark matter. Going to smaller scales, _D$_ grows and eventually overtakes 
the ACDM value due to the increasing contribution of vector perturbations [15]. Thus, the Poisson factor in TeVeS 
shows a characteristic scale dependence and redshift evolution, which should be observable in the correlation ratio 
P(€) defined in section HE. We return to this in section IV D. 

The TeVeS modifications to GR are somewhat more severe than those from the late-time accelerating models f{R) 
and DGP; a simple parametrization in terms of g{k, z) and fciz) does not exist (yet) for this model. Since the TeVeS 
model used here has a slightly different expansion history than ACDM, we compare it with an effective GR-I-DE 
model with the same expansion history. 

D. Comparison with general Dark Energy models 

In the previous sections we pointed out the two main physical quantities affected by gravity, aside from the expansion 

history: the growth factor of matter perturbations, D„i{k, z), and the matter-potential relation, D$_ (fc, z). In modified 
gravity, they are infiuenced by a non-zero metric ratio <?, and a rescaling of the gravitational constant or additional 
degrees of freedom appearing in the Poisson equation. Can these signatures possibly be mimicked by a general Dark 
Energy model ? 

Within GR, anisotropic stress perturbations are the only source of differing cosmological potentials. In order to 
achieve a given metric ratio g{k, z), so that $ + 5' = g^-, one would need components of the dimensionless anisotropic 
stress perturbation H, at late times, of the order of [23]: 

Ilik,z)^^gik,z)^.ik,z). (23) 

Clearly, the components of the Dark Energy anisotropic stress tensor have to be large on small scales in order to 
produce the order unity values of g predicted by modified gravity theories. Any neutrino contribution to H will be 
very small at late times. In addition. Dark Energy density perturbations would also add a contribution of the order 
of Pde^de / {Prn^) to the Poisson ratio V{t} introduced in section HE. While the Dark Energy properties necessary 
to emulate modified gravity seem not to be very natural ones, it is not possible to place stringent constraints on 
general Dark Energy models, in lack of an underlying theory. In the linear regime, the two functions Snsik, z) and 
H(fc, z) are sufficient to emulate any given modified gravity model. However, one might expect that more freedom 
is needed for a Dark Energy model to extend this emulation into the non-linear regime. For example, the f{R) and 
DGP models exhibit a chameleon-like behavior [56] in order to restore GR in high-density environments [37, 39, 42]. 
To emulate this effect, a Dark Energy model would also need a chameleon-like coupling to matter [56]. It would be 
worth studying to what extent general physical constraints can be placed on Dark Energy models which also apply 
in the non-linear regime. 

IV. FORECASTS FOR WEAK LENSING SURVEYS 

We now present quantitative forecasts of the constraining power of future surveys with regard to modified gravity. 
Following the treatment presented in section II, we consider the galaxy-shear and shear-shear correlations, C^*^, and 
C"***, on linear scales {£ < 200 for z > 1). In this regime, the Gaussian error estimate is a good approximation. The 
variance of the cross-correlation C"^^{i), where A, B stand for any of gi, Kj, is then given by: 



riB 



(24) 
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FIG. 2: Left panel: The galaxy-shear cross power C^'^{£)/b for ACDM (black/sohd) and modified gravity theories: f{R) 
(red/dotted), DGP (blue/short-dashed), and TeVeS (green/long-dashed). In case of DGP and TeVeS, the thin lines show the 
corresponding predictions for a GR+DE model with the same expansion history. The foreground galaxies are from bin "F" 
with median redshift of zf = 1.1, and background (sheared) galaxies are from bin "B" with Zb = 3.6 (see text). Right panel: 
Relative deviation of the galaxy-shear cross power C'^"'{£) of modified gravity models from that of a GR-I-DE model with 
identical expansion history, for the same redshift bins as in the left panel. 



Here, /sky is the fraction of sky observed in the survey, nA,B denote the number densities of galaxies per sr in each 
sample, and <7a,b is the shot noise error. For galaxy number counts, ag = 1, while for the shear (convergence), we 
take CTe = 0.35 as the shot noise due to intrinsic ellipticities. We use the exact integrals [equations (7)-(9)] for the 
correlations for £ < 10. For higher £, we use the Limber-approximation expressions [equations (11)-(13)], since the 
deviations are less than 1% and drop quickly for higher £. Note that in any case the bulk of the signal-to- noise comes 
from much higher £. 

In principle one should take into account magnification bias as well, which adds additional fluctuations to the 
observed galaxy overdensity and hence affects C^^ and C^". For a galaxy population with bias b, the observed galaxy 
overdensity dg will be given by [57]: 

6g^b6+{5s-2)K, (25) 

where 6 is the matter overdensity, k is the convergence, and s = dig N/ dm is the logarithmic number count slope as 
a function of magnitude m. Typical values for galaxy surveys are s = 0.2 . . . 0.6 [57, 58]. The galaxy-shear correlation 
will thus be modified as: 

C9''(£) _> c^^i^i) + (5s _ 2) C""^''(£), (26) 

where k/ is the convergence calculated for the foreground galaxy redshift bin. Note that modified gravity will also 
affect the magnification bias term. Since we do not have a realistic estimate for the number count slope in the surveys 
considered here, we do not include this effect. For the galaxy-shear correlation with foreground galaxies centered 
around z — 1.1 and for the range in s given above, we estimate that the effect is around 10% {£ > 50) for ACDM 
and reaches up to 30% for the f{R) cosmology. Note that the parameter s is measurable in surveys. Furthermore, by 
varying the magnitude cut of the galaxy sample, one can vary s to some extent and in this way probe the magnification 
effect. Hence, while straightforward to include for an actual survey, magnification bias will not affect the forecasts 
presented here appreciably. 

Our fiducial flat ACDM cosmology is given by: erg = 0.8, il,,, = 0.27, Ha = 0.73, ilt = 0.046, h = 0.7, = 0.95. 
For the linear matter power spectrum, we use the transfer function from [59]. The parameters used for the modified 
gravity models are stated in section HI. Additionally, we will consider a ACDM cosmology with a non-linear power 
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FIG. 3: Left panel: The shear-shear auto-correlation C'"'{£) for ACDM and modified gravity theories, for high-redshift galaxies 
with median redshift Zb = 3.6. The different colors and line styles correspond to the same modified gravity and GR-I-DE models 
as in Fig. 2. Right panel: Relative deviation of the shear power spectrum C""°(^) from that of GR+DE models with the same 
expansion history, for the same high-redshift galaxies as in the left panel. 



spectrum from halof it [60], in order to assess where the linear approximation breaks down. Note that this is only 
indicative; non-linearities can in principle become important already at larger scales in modified gravity theories, due 
to the non-linear nature of the modified evolution equations for perturbations. 



A. Survey specifications 

We adopt simple survey specifications, as we are mainly interested in the differences between modified gravity 
models from GR in this paper, which do not depend appreciably on the detailed specifics of the survey. The galaxy 
redshift distribution dN/dz is adopted from observations [61, 62] for a magnitude-limited sample with / < 27, as 
parametrized as "Sample I" in [63]. This is roughly what is expected to be attained by LSST [64] or SNAP [65]. The 
median redshift of the distribution is 0.91. We then divide the galaxy distribution into redshift bins defined by z;, Zh'. 

, s „ dN \ ( Z] — z\ ( Zh — zW 

»,W.C-[crfc(-^)-c*(-^)], (27) 

where C is a normalization constant, erfc is the complementary error function, z; and Zh are the lower and upper bin 
boundaries, respectively, and <Jz is the expected photometric redshift error. We adopt = 0.03(1 + z). 

In the following, for studies of the redshift evolution of the galaxy-shear and shear-shear correlations, we divide the 
redshift range from z = — 3 in bins with Az = Zh — zi — 0.4, so that the first bin is defined by zi — 0; z/i — 0.4, 
while the last bin corresponds to zi = 2.4; zu = 2.8. In addition, we define a "background" high-z bin encompassing 
the galaxies from z = 3 to z = 5 (bin "B", with a median redshift zi, = 3.6). We then correlate foreground galaxies or 
shear with the shear of galaxies in bin "B". For studies of the scale- dependence of C^'^{£) and C""'(€), we choose the 
same background galaxies, and use a wider foreground redshift bin defined by z; = 0.8, Zh = 1.6 (yielding a median 
redshift of Zf — 1.1). The constraints on modified gravity do not depend strongly on the number and precise redshift 
of the foreground and background bins chosen. 

We do not attempt to model the galaxy bias in each bin, instead we show correlations divided by 6, where applicable. 
In practice, the bias can be taken as free (scale-independent) parameter to be marginalized ober. We will see that the 
dependence on scale and redshift of the effects of modified gravity should allow them to be disentangled from galaxy 
bias. We will show forecasts for LSST, adopting rig = 50 arcmin"^ as the total observed galaxy density, and /sky = 0.5 
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FIG. 4: Left panel: Relative deviation of the cross power C^'^{£) from that of the corresponding GR+DE models, for the same 
redshift bins as in Fig. 2, but in bins of A£ = 50 with statistical errors expected from LSST. Right panel: Same as the left 
panel, but for the shear-shear correlation of high-z galaxies (as in Fig. 3). 

as the observed sky fraction. Constraints similar to those for LSST are expected from the "wide survey" of SNAP 
[65]. Galaxy-shear and shear-shear correlations offer a large amount of information. As we are mainly interested in 
the specific behavior of weak lensing correlations in the different modified gravity models, we do not attempt to fully 
exploit this information here; instead, we will separately discuss the scale (•^-)dependence and redshift evolution of 
the observables, pointing out generic features of the different model predictions, and whether these effects will be 
observable in future surveys. 



B. Scale dependence 



We first consider the galaxy-shear cross-correlation, correlating foreground galaxies from bin "F" {zf = 1.1) with 
the shear from galaxies in bin "B" (zf, = 3.6). Fig. 2 (left panel) shows C^'^{£)/b for ACDM and the modified gravity 
models, where b is the bias of the foreground galaxies. There are clear differences in the magnitude of the signal, 
with the largest correlation coming from the f{R) model, while DGP and TeVeS yield weaker correlations than the 
corresponding GR-I-DE models with identical H{z) (thin lines). Qualitatively, this is what we expected after the 
discussion of the models in section III (see also Fig. 1). The differences to the predictions of the GR-fDE models with 
the same expansion history are of order 50% or more in the case of f{R) and TeVeS, and of order 10% in the case 
of DGP (right panel of Fig. 2). A scale-dependence of the ratio is also apparent; in the case of the DGP model, the 
deviations become slightly smaller towards smaller scales, while the opposite scale dependence is apparent for the f{R) 
model. This is due to the different values of the metric ratio in the quasi-static (small scale) regime for the different 
models (section III). Note that for DGP and f{R), we expect the GR behavior to be restored on sufhciently small 
(non-linear) scales. In the case of TeVeS, the differences to GR are also sourced by the modified Poisson equation, 
and hence remain present even at small scales. Apparently, for the redshift bins chosen, non-linear evolution becomes 
relevant for £ > 300 in the ACDM case. Hence, any predictions at larger £ should not be taken at face value. We 
found that the deviations from GR in terms of R^'^{£) = C^'^ (t) j ^JCsHl}} ^ which is independent of the linear galaxy 
bias, are very similar to those of C^*^ . 

In Fig. 3 we show the corresponding results for the shear-shear correlation of the high-z "B" bin. The results 
are quite similar to those for the galaxy-shear correlation. The shear-shear correlation offers the advantage of being 
independent of any galaxy bias. However, small-scale modes at low redshifts contribute to the shear-shear correlation, 
so that the effects of non-linear evolution become relevant already at £ ^ 200. In contrast, the contribution to the 
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FIG. 5: Left panel: The galaxy-shear cross power C^'^{£)/b a.t i = 100 for a fixed background galaxy sample with median 
redshift zt = 3.6, as a function of the redshift Zf of foreground galaxies. The different lines correspond to the same modified 
gravity and GR+DE models as in Fig. 2. Right panel: Relative deviation of the galaxy-shear power predicted by the modified 
gravity models from those of the corresponding GR-I-DE models. 

galaxy-shear correlation we considered is concentrated around Zf — 1.1, where the non- linear evolution is somewhat 
less relevant. 

Finally, we determine whether the differences in magnitude and scale-dependence of weak lensing correlations will 
be observable in future surveys. Fig. 4 shows the deviations for C^'^{£) and C'^'^{£) for the same redshift bins as above, 
binned in £ with A£ — 50. Also shown are la statistical errors expected for an LSST-likc survey, using equation (24). 
Apparently, the overall difference in the lensing signal is distinguishable at very high significance, to a comparable 
degree in both galaxy-shear and shear-shear correlations. As both observables are prone to different experimental 
systematics, the possibility of measuring both correlations in the same survey will enable powerful cross-checks of any 
signs of deviations from the GR predictions. Choosing a more conservative background redshift bin from z = 2 to 3, 
with median redshift zi, = 2.4, degrades the expected signal-to-noise of the deviations in C^" and C"" by not more 
than 20%. 

In addition to the overall magnitude, the different scale dependence of C^^{£) and C'^'^{£) in the f{R) and TeVeS 
models should be clearly distinguishable, while the scale dependence of the deviation of the DGP model might be 
difficult to detect. A scale dependence of the deviation of lensing correlations from the GR+DE prediction should 
break degeneracies with other cosmological parameters such as b and erg. The redshift evolution of the lensing signal, 
to which we now turn, can serve as an additional tool for this purpose. 

C. Redshift evolution 

In order to study the redshift evolution of lensing correlations in GR and modified gravity scenarios, we keep the 
background galaxy bin "B" (z = 3 — 5, median Zb = 3.6) fixed while considering foreground redshift bins with various 
median redshifts as explained in section IV A. We keep the width of the foreground redshift bins fixed at Az = 0.4. 
This is only a subset of the cross-correlations possible, which for our purposes serves to show how the evolution of 
correlations differs in modified gravity theories. 

Let us begin with the galaxy-shear correlation, shown in Fig. 5 (left panel) for a fixed £ — 100 as a function of the 
median redshift z/ of the foreground galaxies. The overall behavior as a function of is due to the lensing weight 
function Wl, which peaks at z ^ 0.5, corresponding to half of the distance to the lensed galaxies, of order x(zb). The 
precise evolution with redshift however depends on the evolution of the potentials oc (1 -t- z)Dm{z), which differs in 
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FIG. 6: ie/t panel: The relative difference of the galaxy-shear cross power in modified gravity models, as shown in Fig. 5 
(right panel), but with statistical error bars expected from LSST. Each point corresponds to an independent foreground galaxy 
redshift bin with Az = 0.4; note that errors are correlated as the same background galaxy sample is used for each point. Right 
panel: Same as the left panel, but for the shear power spectrum for a fixed high-z redshift bin with Zb = 3.6 and the same 
foreground redshift bins as in the left panel. 



modified gravity theories (right panel of Fig. 5). We again find the expected result that f{R) shows a larger lensing 
signal, while DGP and TeVeS predict smaUer lensing signals than the corresponding GR+DE models. In all cases, 
we see a significant redshift evolution of the deviations. In the case of f{R) and DGP, the deviations become smaller 
at higher z, as GR becomes restored in the matter-dominated epoch in these models. Thus, ratios of correlations 
such as C^'^{£; Zf = 0.5)/C^'^{£; Zf = 1.5) can serve as probes of gravity. A large part of the observational systematics 
can be expected to drop out in such a ratio, as will the power spectrum normalization erg. However, the different 
galaxy biases of the two foreground samples have to be taken into account. Alternatively, one can measure ratios 
of R^'^ [equation (14)] instead of C^", where similar deviations from GR are seen. It is apparent that at £ = 100, 
non-linearities only become important for very low z foreground galaxies in the galaxy-shear correlation. 

The shear-shear correlation shows a complementary redshift evolution, growing with the foreground galaxy redshfit 
z/, as the amount of mass along the line of sight as well as the geometrical lever arm increase with Zf. The deviations 
in the predicted redshift evolution of the shear-shear correlation for modified gravity models are similar in sign and 
magnitude as those for the galaxy-shear correlation. Again, with its different redshift evolution and set of observational 
systematics, the shear-shear correlation can serve as a cross-check of the galaxy-shear correlation. 

In Fig. 6 we show how precisely the redshift of evolution of C^'^ and C"^"^ is expected to be measurable with a 
survey like LSST. We assume a bin in £ from 50-150, and consider foreground redshift bins spaced by Az — 0.4, so 
that these bins are independent. Note however that the background redshift bin is the same in all cases. It is clear 
that future surveys should be able to resolve the redshift evolution of weak lensing correlations to high precision, and 
distinguish the models discussed here based on this evolution. This also holds when choosing a lower-z background 
redshift bin, e.g. z — 2 — 3, although the redshift range usable for probing the evolution of deviations from GR will 
be restricted to z < 1.5 in this case. 



D. Poisson ratio 



We now discuss forecasts for the measurement of the "Poisson ratio" V{£) = C^'^{£)/C^^ {£) presented in section IIB. 
Fig. 7 shows 'P{£) as a function of £ for the foreground galaxy bin "F" (zf — 1.1) and the background bin "B" (z^ — 3.6). 
It is not completely scale-independent even in the case of ACDM due to the finite width of the foreground galaxy bin 
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FIG. 7: Left panel: The "Poisson ratio" of correlations scaled by the bias, hV{l) [equation (17)], for foreground galaxies with 
mean redshift Zf = 1.1 and (sheared) background galaxies with Zb — 3.6, in ACDM and modified gravity theories. Right panel: 
The same as the left panel, in bins of A£ — 50 and with statistical error bars expected from LSST. 



(Az = 0.4); in section II B we had assumed a (5-distribution in redshift for the foreground galaxies. As expected from 
Fig. 1, TeVeS shows the largest differences in P{£). The modifications to the Poisson equation are very small for the 
assumed f{R) model, while the DGP model shows an offset in 'P{£) due to the different adopted for this model 
(section IIIB). Some of the generalizations of the DGP model considered in [49] have a modified Poisson equation 
which should also leave an observable signature in 1^(1). 

The right panel of Fig. 7 shows the statistical error expected for LSST on the measurement of 'P{£), in bins of 
A£ = 50. For simplicity, the errors on C^^ and C^" where assumed to be independent. In this robust test of the 
Poisson equation, the present TeVeS model should leave an observational signature, as should any modified gravity 
model that shows modifications to the Poisson equation on the order of 5% or greater. Signatures of appreciable Dark 
Energy density perturbations on scales less than ~ Gpc should be detectable in V{i) as well. 



V. CONCLUSIONS 



Uncovering the physics behind the accelerated expansion of the Universe is one of the most compelling open problems 
in astrophysics today. One fundamental question to answer is whether the cause lies in an additional component in 
the energy budget of the Universe, or in a modification of General Relativity on cosmological scales. By considering 
galaxy-shear and shear-shear correlations on large scales, we showed that weak lensing, with its ability to probe the 
scale dependence and redshift evolution of the cosmological gravitational potentials and their relation to matter, can 
serve as a very sensitive probe in discerning between modified gravity and Dark Energy. We focused on the effects of 
modified gravity on the growth of structure, effectively assuming that the expansion history is very well constrained 
through Supernovae la, the CMB, and baryon acoustic oscillation measurements. In practice, distance and growth 
measures should be used jointly to place constraints on modified gravity models. 

Weak lensing is also a sensitive probe of the background expansion history, both through geometry and the growth 
of structure. However, when using the growth of structure to infer the background expansion history (parametrized, 
e.g., by the Dark Energy equation of state w) , one relies on the validity of General Relativity on cosmological scales: 
as we have shown here, modified gravity affects the growth of structure and weak lensing observables independently 
of the expansion history. Fortunately, it is possible to isolate the dependence on the spacetime geometry of lensing 
observables via shear ratios [40, 41]. These techniques will not be affected by modifications to gravity. 

By cross-correlating foreground galaxies with the shear of background galaxies, it is possible to probe the relation 
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FIG. 8; Signal-to-noise (C*^^(^) - C^^+^^(^))/AC"°(^) in £ bins (points) of the deviation of modified gravity (MG) predic- 
tions from GR-I-DE models with the same expansion history, for the galaxy-shear correlation shown in Fig. 4 (left panel) . The 
lines show the corresponding signal to noise of the deviation in the galaxy-CMB cross-correlation via the ISW effect, for the 
same foreground galaxy redshift bin centered at Zf = 1.1. 



between matter and potentials, i.e. the Poisson equation. For this purpose, we introduced the weak lensing correlation 
ratio V{£) [equation (17)], which is, in the hniit of narrow redshift bins, independent of growth effects and the matter 
power spectrum, and thus isolates modifications to the matter-potential relation. However, galaxy-shear correlations 
necessarily depend on the a priori unknown galaxy bias. On linear scales, the degeneracy with the bias can presumably 
be broken by considering the non-trivial scale dependence of the modified gravity signatures; alternatively, one can 
consider the reduced correlation, R^'^{£) [equation (14)]. In contrast, the shear-shear correlation is independent of 
any galaxy bias, while it is more affected by the non- linear gravitational evolution at late times. We stress that 
any modification of gravity should leave a signature in both galaxy-shear and shear-shear correlations, and the two 
methods with their different experimental systematics can serve as cross-checks of the results. 

In this paper, we considered three viable modified gravity models which are, in broad terms, consistent with current 
measurements of the CMB and the expansion history. We compared each model to a GR-I-DE scenario with identical 
expansion history, showing that weak lensing can break this degeneracy even in the case of an extremely tightly 
constrained expansion history. Furthermore, all of these models make definite predictions for weak lensing in the 
linear regime: the f{R) model generically predicts a larger lensing signal than expected in GR; the DGP braneworld 
model predicts a smaller lensing signal than a GR-I-DE model with the same expansion history; and TeVeS predicts 
a weaker lensing signal with considerably modified scale dependence. In all cases, these are robust features of the 
underlying model, linked to the non-zero difference between the cosmological potentials, and modifications to the 
Poisson equation. 

We showed that future wide-field weak lensing surveys, such as LSST and the wide survey of SNAP, can detect 
deviations in weak lensing correlations predicted by these three models with high significance. As an example, Fig. 8 
shows, as a function of £, the signal-to-noise of the deviation of the modified gravity model predictions expected for 
LSST, in the case of the galaxy-shear correlation. For comparison, we also show the corresponding signal-to-noise 
expected for the galaxy-CMB cross-correlation induced by the late-time ISW effect, which is clearly much smaller: 
weak lensing is considerably more powerful in distinguishing modified gravity models from Dark Energy models with 
the same expansion history. As these conclusions hold for the three independent modified gravity scenarios considered 
here, one might expect them to be valid in general for gravity theories that differ significantly from GR on cosmological 
scales. While constraints similar to those of LSST are expected from other proposed surveys, such as the SNAP wide 
survey, these conclusions also hold for less demanding survey specifications. 

Apart from relying on data from future experiments, progress in the theoretical understanding of modified gravity 
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will be crucial in order to improve on the sensitivity of probes of gravity. Using weak lensing to probe gravity on 
smaller scales requires that the non-linear process of structure formation be understood in modified gravity theories. 
Once this is the case, the amount of useful information for probing gravity grows dramatically, say by raising ^max 
from a few hundred to greater than 1000. In addition, one might hope that the degeneracy between modified gravity 
and general, non-smooth Dark Energy models present in the linear theory will be broken in the non-linear regime, at 
least for a class of Dark Energy models restricted by physical constraints on, e.g., the coupling to matter. 

Finally, we point out that the information in weak lensing correlations should be sufficient to place constraints on 
"post-GR" parameters independently of an underlying modified gravity theory, especially for the metric ratio g{k, z) 
and the rcscalcd gravitational constant Gcff . A detailed investigation of this, using all information contained in galaxy 
and shear correlations, is left for future study. 
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